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confined phase. We discuss conceptual and technical problems and present first results of an 
exploratory numerical investigation. 
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1. Introduction 

It is a well known fact, that the Polyakov loop correlator (^P{x)P^ (y) ) can be used to measure 
the static quark-antiquark potential and the correlator of three Polyakov loops tells us about the 
baryon potential. 

As McLerran and Svetitsky showed in [[T]], the expectation value of every combination of 
Polyakov loops is related to the free energy of the appropriate configuration of quarks and anti- 
quarks. Thus, it would be interesting to study the correlation of two Polyakov loops in order to 
determine the diquark potential. 

It has been conjectured that diquarks play an important role in phenomenology, e.g. colour 
superconductivity is motivated by diquark condensates and pentaquarks could be described by 
using diquarks. 

2. The diquark model 

A single quark transforms as a triplet under SU(3) gauge transformations. Two quarks trans- 
form as the tensor product of two triplet states, which can be decomposed as [3] (S> [3] = [3] [6]. 
The antisymmetric [3] -representation transforms in the same way as a single antiquark. This reduc- 
tion of the effective color charge of the two quarks renders the bound state energetically preferable 
and leads to an attractive force. 

There has been a suggestion of how to measure the quark-quark correlation by Jaffe in [Q] that 
has been followed up by other groups [Q]. In case of our study we built up our operators out 
of Polyakov loops only, namly we measured the correlation between two Polyakov loops. This 
method has the advantage that it is not necessary to fix the gauge. Furthermore we do not need a 
third quark line like in Jaffe's suggestion. For alternative approaches to measure the free energy 
see also [|], [g] and ||]. 

3. Diquark free energy 

On the lattice, the Polyakov loop expresses the propagation of a single heavy quark along a 

closed loop in periodic imaginary time: P{x) =Tr I Yl Ut{x + ki) \ . 

\k=i ) 

McLerran and Svetitsky showed [[T]] that the free energy of a configuration of n quarks and n 
antiquarks at positions ;ci , . . . ,;c„ and ji , . . . ,^77 respectively is related to the expectation value of the 
appropriate combination of Polyakov loops as 

(P(X1 ) • • • PiXn )P^ (jl ) • • • ijn) ) = exp {-aN, {Fnqflq " Fqo)] • O- 1) 

Here, a is the lattice spacing, Nt is the lattice extension in the time direction and Fqo is the free 
energy of the vacuum without any quarks. Especially for a quark-antiquark pair, this equation 
leads to the static quark-antiquark potential which is the free energy at zero temperature 

{P{x)P\y)) = tx^\-aN,{Fq-,{\x-y\)-Foo)] . 
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From equation (3.1), one can construct the normalized free energy of a specific configuration, 
i.e. the gain in the free energy of a finite configuration of quarks relative to a configuration where 
all the quarks ai^e infinitely separated. 
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Especially for two quarks, this relation leads to the diquark free energy 

The single quark free energy of the vacuum Fq — Fqq is included into the operator in order to remove 
the self energy of the two quarks. 



4. The Polyakov loop in pure gauge theory 

It is well known, that a pure SU(3) gauge theory ex- 
hibits an additional Z(3) center symmetry. The action is 
not changed if all the time links connecting two neighbor- 
ing timeslices are multiplied by the same element of Z(3). 
The Polyakov loop however does change under such a 
transformation 



P(x) 



■P{x) {n = 0,1,2). 
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Figure 1: quenched Polyakov loop density 
distribution from a 6^ lattice in the confine- 
ment region 



This makes the Polyakov loop sensitive to the cen- 
ter symmetry. Thus, in the quenched theory, it serves 
as an order parameter for the spontaneous breakdown 
of this symmetry, which is associated with the confine- 
ment/deconfinement phase transition. 

Since the update algorithm will generate configura- 
tions from all three Z(3) sectors with equal probability, 

the values are distributed symmetrically in all of the three sectors as can be seen in figure |. Thus 
the loops of the different Z(3) sectors will cancel out each other, which results in an exactly van- 
ishing expectation value for the Polyakov loop. Therefore, in order to measure a non-vanishing 
Polyakov loop, one has to break the center symmetry. 

There have been some propositions how to break this symmetry. For example, one can take 
|P(x)|, iyRe[P(x)^], rotate all values around zero so that the phase of each loop lies in the range 
from — J to J like in or just ignore all loops which do not lie in the sector around the positive 
real axis as has been done in |j8|]. But all these methods seem somehow ai^biti'ary. This differs from 
the case when dynamical fermions ai^e included into the simulation. In the unquenched theory, the 
fermion determinant exphcitly breaks the center symmetry. 
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5. Simulation details 

In full QCD, the Polyakov loop does not serve as an order parameter any more, but equation 
( p7l| ) still remains valid. Because of the exponential decay, it is an expensive task to determine 
the Polyakov loop and its correlators accurately for large lattice extension in the time direction or 
equivalently for low temperature. But the low temperature region is especially interesting. One 
drawback of the dynamical staggered fermions in full QCD is the non-applicability of exponential 
error reduction such as the Luscher-Weisz multilevel algorithm [^. Therefore, the lattices which 
can be studied are restricted to small time extension. 

We were able to study lattices with a time extension up to A'j = 6 and a spatial extension of 
Ns = 18. The first calculations were done on a 4^ and a 6^ lattice with the Wilson plaquette action 
and 4 flavours of staggered fermions with a mass of niqa = 0.05 in lattice units. 

With this data set, we calculated the di- 
quark free energy at different distances as a 
function of /3 . The result for r = a is shown 
in figure As can be seen from this figure, 
there is a clear correlation between the two 
Polyakov loops. The sharp jump at j3 5 
corresponds to the confinement/deconfine- 
ment phase transition. 

In order to determine the quark-quark 
correlator for larger distances, we did the 
same calculations on a 6^^ lattice at j8 = 5. 14 
which is still in the confinement region and 
corresponds to a lattice spacing of a = 0.28 
fm. The form of the qq free energy in com- 
parison to the quark-antiquark potential is 
shown in figure ^ One can see, that the free 
energy falls off at shorter distances which 
indicates an attractive force. 

It is preferable to determine the form 
of the qq free energy for larger separations. 
Since the Polyakov loop decays exponen- 
tially with Nf, as described by equation (^), we kept the temporal lattice extension fixed to A'; = 6 
and increased the spatial lattice size. 

The following results were obtained from a second data set, which was generated using 2-1-1 
staggered fermions with stout smeared links [|l^] and physical quark masses and the Symanzik tree 
level improved gauge action on an 18^ x 6 lattice. We used the finite temperature configurations 
which were generated in 

Figure ^ shows the diquark free energy as a function of the distance between the Polyakov 
loops and compare it with the quark-antiquai^k free energy. Note that our data sets correspond 
to finite temperature samples with T = 0.967^ and T = 1.14rc. Similar to the quark-antiquark 
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Figure 2: qq free energy from a 4^ lattice 
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Figure 3: qq and qq free energy from a 6^ lattice at /3 
5.14, the errors are smaller than the dots. 
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Figure 4: quark-quark free energy compared with the quark-antiquark free energy for ^ = 0.96 (left) and 
Y = 1.14 (right). In the left plot, a box-averaging has been used for r > 0.5 fm to reduce the errors. The 
boxes were smaller than 0.04 fm. In the right plot, the errors are smaller than the dots. 



attraction, there is an attraction between two quarks. The quark-antiquark attraction is significantly 
stronger than the quark-quark attraction. For our temperatures and for r > 0.8 fm both the qq and 
the qq signal are compatible with zero. 

In agreement with the results of a calculation in the quenched theory presented in [ |l2| ] [0], we 
see a flattening of the qq free energy above Tc which becomes constant at about the same distance. 
We also see a flattening below Tc, whereas in the quenched calculations, neither the qq nor the qq 
free energy seem to approach a constant value below Tc. 



6. Conclusions 



We see a clear signal for a Polyakov-Polyakov loop correlation in full QCD without any need 
to fix the gauge. The diquark free energy has a similar form as the quark-antiquark free energy. 
This is a sign for an attractive force between two quarks, but it is weaker than the attraction between 
a quark and an antiquark. 

Accurately calculating the Polyakov loop expectation value is still an expensive task. In par- 
ticular, the investigation of zero temperature lattices is restricted by the exponential decay of the 
Polyakov loop. As a consequence, it is not yet clear, how this Polyakov-Polyakov coiTclation can 
be determined in the limit T — > 0, which is of physical interest. 
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